Quantum theory of excess noise 
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We analyze the excess noise in the framework of the conventional quantum theory of laser- 
like systems. Our calculation is conceptually simple and our result also shows a correction to the 
semi-classical result derived earlier. 
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In attenuators, amplifiers and laser systems, the noise of the signal is increased by the interaction with a reservoir. 
In usual laser systems this is reflected in the formula for the Schawlow-Townes linewidth which provides a simple 
relation between the linewidth and the gain or loss of the system. Under certain conditions, however, the noise entering 
from the reservoir can exceed this minimum amount by a large factor — the so-called excess-noise or Petermann factor. 
The appearance of this factor was predicted by Petermann [Q in the context of gain-guided semiconductor lasers. 
Within a semi-classical theory |^| , this concept was later generalized to other systems in quantum electronics. 

However, for a general system a complete quantum mechanical derivation is still lacking. Most of the considerations 
■ so far include the quantum noise properties ad hoc into an otherwise classical theory. Only a few simple systems have 
\ been discussed quantum mechanically in a rigorous way |]Q . 

During recent years, many interesting experiments have been carried out: The first experiments used cavities with 
large output coupling || enhancing the noise by a few times only. However, geometrically unstable laser cavities show 
excess-noise factors up to a few hundreds for both solid state lasers |)| and gas lasers [|o| . Also the coupling of the 
polarizations of the laser light [[ll] and the insertion of a small aperture |l2| can lead to large excess noise. 

We derive the noise properties of an amplifier starting from the usual quantum mechanical description of the 
radiation field. First we derive a multi-mode master equation for the amplification. This is done in a way analogous 
to the maser theory assuming a reservoir of excited atoms Q. Then we are able to define quasi modes of the total 
system which show noise enhanced by the so-called isT-factor. Our calculation is conceptually simple and transparent; 
it applies to all linear systems. We derive a correction to the semi-classical theory of Ref. ||. We subsequently 
generalize our result by including damping which allows us to apply it to the case of a laser below or up to threshold. 
We explain the main features of the excess noise and discuss its properties. For the experimentally relevant cases, we 
show that our calculations justify the use of the ordinary semi-classical treatment of excess noise. Thus we refrain from 
providing any explicit physical realization of our abstract model. Such a calculation would, in all essential features, 
reproduce the calculations carried out in the various model systems discussed in the literature ||. In this paper we 
O^ 1 do not include the effects of nonlinear saturation which needs to be considered separately 

Here we briefly recall the standard expressions of the quantized electromagnetic field in order to define our notation. 
We use orthonormal real mode functions u n (x) of the electromagnetic field with frequency uj n which fulfill the boundary 
conditions for the given configuration in the whole "universe" and satisfy the orthonormality relation 
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d xu n (x)u m (x) = S nm , (1) 



where V is the volume of the whole space. Note that the mode function u n (x) is a vector including the polarization 
orientation and that we choose them to be real for convenience. The electric field operator [[l3| then reads 

E{x) = ^ £ nU n (x) (a n + a*) , (2) 

n 

where a n and a) n are the usual creation and annihilation operators of the field excitations and the so-called vacuum 
field amplitude is 



(3) 

We assume an amplifier medium consisting of two level atoms. This model has frequently been used for the quantum 
treatment of laser systems . The atoms start in the upper level, interact independently with the field for a short 
time before they are repumped incoherently to the upper level. The master equation for the density operator of the 
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field is then obtained by tracing out the atomic variables. To describe the interaction between the field and a single 
atom at position x, we use the interaction Hamiltonian 



H = rf^ e n u n (x){a) n a + a„<7 + ), (4) 

n 

where are the raising and lowering operators for the two atomic levels coupled by the dipole moment d. After a 
short interaction time r the change of the reduced field density operator p(t), in the interaction picture, is given by 



t 2 r , -i 

Sp(t) = ^2 X! £ n£ m [u n {x)d][u m (x)d]e l ^ n ^ ulm)t p],p(f)a ra - a m a\p{t) - ptyama^j 

+0(r 3 ) (5) 



for one atom. We now assume that the atoms are introduced uniformly distributed over a volume V ^ V in the 
upper level at a rate R. In principle we could let R depend on the position x and the reservoir average could include 
an average over the orientation of the dipole moment d. Without such refinements, the rate of change of the ensemble 
averaged density operator is given by 

p(t)=R±; J d 3 x5Mt) 

= 2 £ I ™,« e,K " Um)i {2aif>(t)a m - a m alf}{t) - p(t)a m alj , (6) 

with the matrix elements 

£m,n = -^-£n£m-^ J d 3 x [u n (x)d}[u m (x)d]. (7) 
Transforming from the interaction picture to the laboratory frame we get 



dt f>{t) = i^^L m . n {2a\ l p(t)d m -a m a) n p(t) - p(t)o m a^} - i w n [al a n , p(t)}. (8) 



d 



n.m 



This multi-mode master equation with a symmetric matrix L m>n allows us to define the quasi modes and derive their 
noise properties in the following. Its important feature is the coupling between modes with different frequencies due 
to the reservoir. 

We may use the same kind of model also to derive dissipative losses of the field if we assume a fraction of the 
reservoir atoms to be initially in the lower state. These can be introduced into the system with a distribution differing 
from that of the amplifying atoms. But we will focus now only on the linear amplifier and discuss the case with 
damping later. 

We are looking for a mode operator A which obeys 

±{A} = (±-iQ){A) (9) 

for an arbitrary field state. Here f2 is the real frequency and A is the real amplification rate. We write this mode 
operator in terms of the free field mode operators as 

£A^^e n c n d n (10) 

n 

with the expansion coefficients c„. This transformation includes the vacuum-field amplitudes e n and we define 



£ = J 2^7: because then the classical field amplitudes e n (d n ) obey the same transformation. Inserting Eq. (^|) into 
(|j"o|) we get an eigenvalue equation 



^ ' ( 2^ m n — i^n,m^n J ~ — c n — iQ)c m (H) 
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for the non-Hermitian matrix L mj „ = (^I m „ — iS n)m (jj n )^. Here Cn^ is the right eigenvector of L m . n \ the corre- 
sponding left eigenvector is e^Cn @]- The superscript v distinguishes the different eigenvectors. 

The only properties of the left and right eigenvectors of non-Hermitian matrices which we need for our analysis are 
their mutual orthogonality and completeness : The eigenvectors fulfill the orthogonality condition 

^cl^^E^ 2 (12) 

n n 

and the completeness relation 

_ / j>)>) \ 

E " ,2 = ^ m (13) 

( ) 2 

with ^2 n ' £ n' c n r 7^ 0. We can now uniquely define the set of quasi-mode operators as 

n 

with the vacuum field amplitude 



The inverse transformation is 



£ " = ^- (15) 



(16) 

" Em £ m c ™ 



Consequently the positive frequency part of the electric field operator is given by 



2 Jy) 



The quasi-mode eigenfunctions 



satisfy an orthogonality relation 



= Y J £vU„{x)A„. (17) 



Uu{x)=Y< — t ^ n{x) (18) 
71 Z^m t m L -m 
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d°xU v (x)U M (x) = 8 Vtll (19) 
with their adjoint quasi-mode functions 

U v {x) =J2 c n ) Mx)- (20) 



The quasi-mode functions have the norm 



nI = i / d 3 z [/,(*){/; (x) = E " 4|c ^ ) |2 

2—ira c m c " 1 
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(21) 
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and their adjoints 



»|2 



(22) 



From now on we only consider one quasi mode and drop the index v. This is well justified, because, after a long 
enough time, that quasi mode which has the largest amplification rate A is dominating in the sum Eq. (17). Thus all 
expectation values can be calculated with just this largest contribution. 

Later we need the properties 



n 



En £ n\ c n\ 2 



En £ nl c "| 2 



and 



(23) 



(24) 



which can be obtained from the real and imaginary parts of Eq. ( |ll| ) after taking the scalar product with the vector 

2 „* 



e^Cm- Note that Q is the mean frequency with respect to the probabilities 



Pn 



4M 2 



(25) 



En' £ n' l C ™' 1 2 

We calculate now the noise of the slowly-varying quadrature operator 

X(x) = £(U(x)Ae iQt + U*(x)A^e- int ) (26) 
in a frame rotating with Q. After some straightforward calculation we get the time evolution of the noise 

d 



dt 



(AX(x)f = X(AX(x)) 2 + \\U(x)\ 2 J2 £ lK\ 2 . 



We have to compare this with the noise of the usual single-mode-amplifier master equation 

^p[t) = ^ {2a)p(t)a - aa)p{t) - p(t)aa f } - ifl[a f a, p(t)} 



(27) 



(28) 



for the mode frequency f2, normalized mode function u(x), and the vacuum field amplitude £ — J 37^7- I n this case 
the noise is given by the equation 



2 c 2 



dt 



{AX(x)Y = X(AX(x)Y + \\u(x)\ 4 £ 



(29) 



Hence, after averaging over position in Eqs. fl27| ) and d29|) and using Eq. (|2l|), we can read off the excess- noise factor 

K 



|2_E4|cJ 2 E„4|Cn, 

2 / , £ nl C " 



£ 



£ 2 lEm £ m C m 2 \ 



T, £ lK\ 



Em £ rn c n 



(30) 



From the triangular inequality V]£^|c n | 2 > IEm £ m c m 2 | follows K > 1. The peculiar mode coupling of the non- 
Hcrmitian eigenvalue equation (11) causes this enhancement of the reservoir noise entering the quasi mode. 
We can now compare the result, Eq. (|30"|), with the expression 



K 



J d 3 xU(x)U*(x) J d 3 xU{x)U*(x) 



N N = 



Jd 3 xU(x)U(x)\ 

Era l C "| 2 En £ n\ c n\ 2 
I Em £ m c m 2 \ 



(31) 
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which has been derived from a semi-classical theory [fjj. We see that they are almost identical when we use Eqs. (JTs|) , 
( pl| ) and (|2^), The ratio of the factors is given by 



K _£ 2 Zn\Cn\ 2 _EJCn| 2 £„k| 2 £ 4 



# £ n 4M 2 (EJc„l 2 4) 2 



n(l/u)wl + 0(Aw/O)' (32) 



with the mean inverse frequency lu^ 1 with respect to the probabilities p n introduced above in Eq. (|25|). Here we have 
used the definitions of e„, £ and fl in Eqs. (|J), ([l5]) and (^3|). The quantity Aui is a measure of the bandwidth of the 
quasi mode in terms of the modes of the universe. From Eqs. © and (|l]) we see, by using Schwarz' inequality, that 
K < K. The correction to the semi-classical result is small for the optical frequency domain where the bandwidth 
(< 10 10 Hz) is negligible with respect to the mean frequency (> 10 14 Hz). But in the micro-wave regime (GHz) the 
correction may be essential. 

We note that if the quasi-modc frequency f2 was defined with respect to the weights |c„| 2 instead of Eq. (25), the 
ratio Eq. ( |3^ ) would be unity. Indeed, |c„| 2 is the weight of the classical field amplitudes e n (a n ) in Eq. (10) and 
it is the weight of the universe modes u n {x) in the adjoint quasi- mode functions U u (x) in Eq. (|20|). On the other 
hand, the weight of the universe modes in the quasi-mode functions in Eq. ( |l8| ) is £^|c„| 2 . But, Eq. ( p3[ ) shows that 
the geometric mean of the two classical possibilities gives the proper distribution p n for the frequencies with the 
eigenfrequency Q as mean value. 

As a next step we could derive damping for the quasi modes by considering the limit of a continuum of modes of the 
universe. In the case of classical field modes, this was done by Lang, Scully and Lamb fl^| in one dimension for a cavity 
with one perfect and one semi-transparent mirror. It is interesting to note that in this special case, the authors could 
explicitly prove that only one quasi mode exhibits amplification whereas all other quasi modes experience attenuation. 
However, in the context of discussing the properties of excess noise in this paper, we are satisfied by adding damping 
due to a separate reservoir. In an analogous way as before we can derive the multi-mode master equation 

= \ Lrn - n {^n/^K" ~ a m alp(t) - p(t)a m al} 



n.m 



7jE r " l >" {2a n p{t)a\ n -a] n a n p(t) - p(t)a] n a n } - i^2u n [ala n , p(t)}. (33) 



n , m 



The symmetric matrix r mj „ is defined similar to the matrix L m ^ n in Eq. (0), with a possibly different volume V and 
a possibly different atom injection rate R for the lower state. 

We are interested in the case when the strength of the amplification is adjusted such that the quasi-mode amplitude 
obeys the equation 

!<i) = -ffi(i) (34) 

yielding a stationary modulus. At threshold, the laser quasi-mode is oscillating without change of amplitude. Thus, 
the effects of amplification and attenuation must compensate exactly on the average, but only for this quasi-mode. 
This does not mean that the two contributions cancel at the level of the master equation, but only that we can find 
an eigenmode with a purely imaginary eigenvalue. In fact, assuming exact cancellation in the master equation leads 
us to a trivial case, see below. In Ref. || it is assumed that the same applies to the saturated gain at steady state 
operation above threshold; here we consider only the linear behavior close to threshold. We return to the discussion 
of the nonlinear problem in a forthcoming publication. 
These assumptions lead to the eigenvalue problem 



E 

n 



We thus require the amplification rate A as given in Eq. (gj) and the damping rate 

7 = E„4|c n | 2 (36) 

to be equal; 7 = A. Proceeding as before and averaging over position, we obtain the equation 
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-^(AA) 2 = (7+ \)£ 2 K = 2\£ 2 K 



(37) 



for the noise with the same excess- noise factor K as in Eq. (p0|) . This describes a diffusion process with the diffusion 
constant 2Dx = 2\£ 2 K. The lincwidth of the laser at threshold |j] [3| is then given by the phase diffusion constant 
2£> = 2D X /I = K\ 2 /AP with field intensity I = 4£ 2 A 2 (iU) and output power P = 7//(4£ 2 ) = 7 A 2 (it A) 
measured in photon energies. The measurement of the laser intensity is not straightforward in the case of a quasi- 
mode. However, if the detector is mode matched to the outgoing quasi-mode profile, the maximum intensity I is 
obtained. All other detector arrangements will miss some intensity and give a lower value. In the present case, we 
find the noise to be purely due to phase diffusion. This derives from our assumption that we consider the steady state 
near threshold, not a linearized noise theory around the saturated gain. In the latter case, we additionally would 
expect a line width contribution from the amplitude fluctuations. 

We consider now the steady state case when the mode coupling is mainly caused by the damping, which holds if 
L m ,n ~ ^m.n- Then we get from Eq. (|35| ) the eigenvalue equation 




which is the equation describing a system without an amplifying medium. Since with this assumption the eigenvectors 
are the same for Eqs. ( |35|) and ([38]), also the if- factor is the same in both cases. This result justifies the commonly 
used procedure to calculate the if-factor of a laser system from the cavity-decay properties alone. Our assumption 
probably gives an over-simplified picture of a real laser. However, in most systems the gain medium is distributed 
uniformly over an essential part of the mode volume. From the definition, Eq. (0), then follows that L m n ps A5 m ,„. 
Besides this, the gain is often due to a resonant interaction whereas the damping is highly broad band. In that 
situation it seems reasonable to assume the gain to be diagonal in the "universe" modes and the mode-mode coupling 
to derive mainly from the loss mechanisms. These contain diffractive losses at the cavity edges and transmission 
losses, and hence their distribution is bound to differ greatly from the distribution of the amplification. In comparison 
to the off-diagonal elements of r m .„, the off-diagonal elements of L m n can be neglected. 

From Eq. (^) we can also see in which case the minimum excess noise K = 1 can be achieved. For L m>n = r m „ we 
have a n — A n . Losses and amplification are compensated in each volume element. This means that, when damping 
and amplification are acting equally in the same volume V, the true modes and the quasi modes are identical. This 
can happen when losses are mainly provided by the same entities which give the gain like in dye lasers. In most 
lasers, however, it is mainly the smallness of the off-diagonal elements of the matrices L m%n and r m>n which leads to 
a A-factor close to one. 

Summarizing our main results, we have derived the excess-noise factor K within the framework of the conventional 
quantum theory of laser-like systems. Our calculation is conceptually simple and our result also shows a small 
correction to the well established semi-classical one. 
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